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We study the evolution of matter density perturbations in Galileon cosmology where the late- 
time cosmic acceleration can be realized by a field kinetic energy. We obtain full perturbation 
equations at linear order in the presence of five covariant Lagrangians Li (i = 1, ■ ■ ■ , 5) satisfying the 
Galileon symmetry d^tf) ~~ > d^cf) + 6 M in the flat space-time. The equations for a matter perturbation 
as well as an effective gravitational potential are derived under a quasistatic approximation on 
subhorizon scales. This approximation can reproduce full numerical solutions with high accuracy 
' for the wavelengths relevant to large-scale structures. For the model parameters constrained by 

£\j i the background expansion history of the Universe the growth rate of matter perturbations is larger 

than that in the A-cold dark matter model, with the growth index 7 today typically smaller than 
0.4. We also find that, even on very large scales associated with the integrated-Sachs- Wolfe effect in 
cosmic microwave background temperature anisotropies, the effective gravitational potential exhibits 
a temporal growth during the transition from the matter era to the epoch of cosmic acceleration. 
CNJ ■ These properties are useful to distinguish the Galileon model from the A-cold dark matter in future 

<N ' high -precision observations. 



I. INTRODUCTION 



o 

The large-distance modification of gravity has received much attention as a possible explanation for the cosmic 
£Li acceleration today [l|. Many modified gravitational models of dark energy have been already proposed-including 
q those based on f(R) gravity @, scalar-tensor theories Q, Gauss-Bonnet gravity and its generalizations Dvali- 
$_i ■ Gabadazde-Porrati (DGP) braneworld Q, and Galileon gravity |(|. In f(R) gravity, for example, the viable models 
' Q are designed to have a large mass of a scalar gravitational degree of freedom ("scalaron" Q) in the regions of 
high density for the compatibility with local gravity experiments. In this case, as long as the so-called chameleon 
mechanism |9|is at work, the interaction between the scalaron and baryons can be suppressed to satisfy local gravity 
, constraints [10| . 

J> 1 There is another mechanism to decouple the fifth force from baryons at short distances even in the absence of the 
! scalar- field potential. Nonlinear effects of field self- interactions can allow the recovery of General Relativity (GR) 
inside a so-called Vainshtein radius [TlJ. In the DGP model a longitudinal graviton (i.e. a brane-bending mode </>) 
gives rise to the self-interaction of the form (r^/m p i) □0(9 A1 <^<9 A1 </>) through the mixing with a transverse graviton, 
where r c is a crossover scale of the order of the Hubble radius Hq 1 today and m p i is the Planck mass [l2|. In the local 



region where the energy density p is much larger than r c 2 m^ l the nonlinear self- interaction can lead to the decoupling 
of the field from matter. However the DGP model suffers from a ghost problem in addition to the difficulty for 
consistency with the combined data analysis of Supernovae la (SN la) and Baryon Acoustic Oscillations (BAO) fl4j . 

The self-interacting Lagrangian □0(d A ' (f>d '^(fi) appearing in the DGP model satisfies the Galileon symmetry d^ifi — > 
dfj,(j) + b^ in the Minkowski background. Imposing the Galileon symmetry in the flat space-time one can show that the 
field Lagrangian consists of five terms C\, ■ ■ ■ , £5, where the term □(/>(<9 A '</>9 M 0) corresponds to £3 Q. In Refs. [l5[ 
5_i ' these terms were extended to covariant forms in the curved space-time. Moreover one can keep the equations of 
5^ ! motion up to the second-order, while recovering the Galileon Lagrangian in the Minkowski space-time. This property 
is welcome to avoid the appearance of an extra degree of freedom associated with ghosts. 

In Refs. [HI, [TtJ two of the present authors studied the cosmological dynamics of covariant Galileon theory in the 
presence of the terms up to £5 (see Refs. [l8|]-[3!| for related works). There exist de Sitter (dS) solutions responsible 
for dark energy driven by the field kinetic energy. Refs. (Tfjl . also clarified the viable model parameter space in 
which the appearance of ghosts and instabilities of scalar and tensor perturbations can be avoided. 

In the covariant Galileon cosmology the solutions finally converge to a common trajectory (tracker), which is 
characterized by the evolution </> oc H^ 1 (H is the Hubble parameter) [l6|. The epoch at which the solutions 
approach the tracker depends on the initial conditions of the variable r\ = (j) ( isHds/(<i ) H), where the subscript "dS" 
represents the values at the dS point. For smaller initial values of T\ the approach to the tracker, characterized by 
r\ = 1, occurs later. In Ref. [40] it was shown that the combined data analysis of SN la, BAO, and the cosmic 
microwave background (CMB) shift parameter tends to favor a late-time tracking behavior around the present epoch. 
This comes from the fact that the dark energy equation of state wde — — 2, which corresponds to the one for the 
tracker during the deep matter era [3, is difficult to be compatible with a number of observations. 
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In order to constrain the Galilcon model further, it is important to know the evolution of cosmological perturbations 
from the matter era to the epoch of cosmic acceleration. In particular, the modified growth of matter perturbations 
S m affects the matter power spectrum as well as the weak lensing spectrum [4l| . Moreover the modification of gravity 
manifests itself for the evolution of the effective gravitational potential $ c ff related with the integrated-Sachs- Wolfe 
(ISW) effect in CMB anisotropics. In this paper we shall derive the equations of cosmological perturbations in the 
presence of the five Galileon terms <C, (i = 1, • • • ,5) and numerically integrate them to find observational signatures 
of the Galilcon model. We also obtain convenient forms of the effective gravitational coupling G e g and the anisotropic 
parameter r\ between two gravitational potentials, under a quasistatic approximation for the modes deep inside the 
Hubble radius. 

This paper is organized as follows. In Sec. [II] we review the Galileon cosmology in the flat Friedmami-Lcmaitre- 
Robertson- Walker (FLRW) background. In Sec. IIIII the equations of cosmological perturbations in the Galileon 
model are derived in the presence of nonrelativistic matter. In Sec. IIVI we obtain the equations for 6 m and $ c ff 
under a quasistatic approximation on subhorizon scales. In Sec. |V]we present numerical results for the evolution of 
perturbations by integrating the full equations of motion for the wave numbers relevant to large-scale structures and 
CMB anisotropics. Sec. IVIl is devoted to conclusions. 



II. BACKGROUND GALILEON COSMOLOGY 



Let us consider the following action 

- M 2 1 5 



/M l 
d 4 xy/^g~ ~2~ R + o 51 ClCi + Pm s 
i=i 



(1) 



where g is a determinant of the space-time metric g^, M p \ = (87TG)" 1 / 2 is the reduced Planck mass (with G being 
gravitational constant), R is a Ricci scalar, and Ci are constants. The five covariant Lagrangians Ci (i = 1, ■ ■ ■ ,5) are 
given by [lj| 

A=M 3 ^, £ 2 = (V</>) 2 , C 3 = (□<£) (V</>) 2 /A/ 3 , 
U = (V^) 2 [2(Q/>) 2 - 2<t>,^ v - i?(V0) 2 /2] /M 6 , 

C 5 = (V^) 2 [(D</>) 3 - 3(D0) <t>;^ v + 2^ v - H-^^GvpMM* , (2) 

where M is a constant having a dimension of mass, and G vp is the Einstein tensor. These Lagrangians respect the 
Galilcon symmetry in the Minkowski space-time. Moreover the field equations are kept up to the second-order in time 
derivatives. The last term, p m , represents the pressure of a perfect fluid whose energy density and equation of state 
are p m and w = p m /pm, respectively. The pressure depends on a chemical potential /i = (p m +p m )/n (n is a number 
density) and an entropy per particle s. 

Since we are interested in the evolution of matter density perturbations long after the radiation-domination, we 
take into account a nonrelativistic fluid (w ~ 0) only in the following discussion. In the FLRW space-time with the 
scale factor a(t), the variation of the action @ leads to the following background equations fl6l Il7|: 

3M^H 2 = PDE + Pm , (3) 
3M^H 2 + 2M* 1 H = - PBE , (4) 

where H — a/ a is the Hubble parameter (a dot represents a derivative with respect to cosmic time t), and 

Pde = -cii\/ 3 0/2 - C2</> 2 /2 + 3c 3 H<j) 3 /M 3 - 45c4H 2 <p 4 /(2M 6 ) + 2lc b H 3 ft / M 9 , (5) 
p DE = Cl M 3 <j)/2 - c 2 ^ 2 /2 - c 3 (j> 2 4>/M 3 + 3c 4 ft[&H4> + (3H 2 + 2H)(j)]/(2M 6 ) - 3c 5 ii^ 4 [5iJ0 + 2(iJ 2 + H)^]/M 9 .(6) 

The nonrelativistic fluid satisfies the continuity equation p m + 3Hp m = 0. 

We shall consider the case in which the late-time cosmic acceleration is realized by the field kinetic terms, i.e. 
Ci = 0. There is a de Sitter (dS) solution characterized by H = H^s = constant and <fi = (f>ds = constant. From 
Eqs. © and (g]) it then follows that 

c 2 x 2 ds = 6 + 9a - 12/3 , (7) 
c 3 x 3 dS = 2 + 9a - 9/3 , (8) 
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where x ds = 0ds/(-HdS-Mpi) and 

a ee c 4 xi s , /3 ee c 5 x 5 ds . (9) 

We normalize the mass M to be M 3 = M p iH^ s . Since -ffdS is of the order of the present Hubble parameter (H « 
10- 60 M pl ), we have that M « 10 _40 Mpi. 

In order to discuss the cosmological dynamics, it is convenient to introduce the following variables 

_ 4>dsH dS _ 1 ( <p \ 

n = — r— — , r 2 = — ^— . (10) 



<f>H n ^ dS y 

At the dS point we have r\ = 1 and r 2 = 1. The Friedmann equation can be written as £^de + ^ m = 1, where 
n m ee p m /(3M^H 2 ) and 

" DE = 3Ji/°ff2 = ~5 (2 + 3a ~ 4/3)r ' r2 + ( 2 + 9a - - yonra + 7/3r 2 . (11) 

We also obtain the autonomous equations 

r[ = i(n-l)ri[r 1 (ri(-3a + 4/3-2)+6a-5^)-5 / 8] 

x [18 + 3r 2 (r?(-3a + 4/3 - 2) + 2rf (9a - 9/3 + 2) - 15na + 14/3)] , (12) 

r 2 = -^-[r 2 (6r 2 (r 2 (45a 2 - 4(9a + 2)/3 + 36/3 2 ) + 7(9a - 9/3 + 2)) + r?(-66(3a -4/3 + 2) 

-3r 2 (-2(201a + 89)/? + 15a(9a + 2) + 356/3 2 )) - 3ria(123r 2; 3 + 36) + 10/3(21r 2 /3 - 3) 

+3?iV 2 (9a 2 - 30a(4/3 + 1) + 2(2 - 9/3) 2 ) + 3r?r 2 (3a -4/3+ 2) 2 + 3rfr 2 (9a -9/3 + 2) (3a - 4/3 + 2))], (13) 

where a prime represents a derivative with respect to N = In a, and 

A ee 2r^r 2 [72a 2 + 30a(l - 5/3) + (2 - 9/3) 2 ] + 4r 2 [9r 2 (5a 2 + 9a^ + (2 - 9^)/3) + 2(9a -9(3 + 2)] 

+4r?[-3r 2 (-2(15a + 1)0 + 3a(9a + 2) + 4/3 2 ) - 3a + 4/3 - 2] - 24ria(16r 2/ 3 + 3) + 10/3(21r 2 /3 + 8) . (14) 

The Hubble parameter obeys the following equation 

* = _*A_iL. (15) 

# 4n 4r 2 V y 
We define the dark energy equation of state wde and the effective equation of state w e g, as 

»DE= , W eS = -1 - — — . (16) 

/ODE oil 

Using the continuity equation pue + 3H(pnE + .Pde) = 0, wc obtain the relation wbe = WeS — ^de/(3^de)- 

There are three distinct fixed points: (A) (ri,r 2 ) = (0,0), (B) (ri,r 2 ) = (1,0), and (C) (ri,r 2 ) = (1, 1). As we see 

from the definition in Eq. (|TU| the point (C) corresponds to the dS solution, which is always classically stable fl7| . 

The point (B) is a tracker solution found in Ref. [l6j], along which the field velocity evolves as <fi oc 1/H. The fixed 

point (B) is followed by the stable dS solution once r 2 grows to the order of 1. Depending on the initial conditions of 

ri, the epoch at which the solutions approach the tracker is different. 

In the following we summarize the background evolution in two regimes: (i) r\ -C 1, r 2 -C 1 and (ii) r\ = 1 [III [13 ■ 

• (i) n < l,r 2 < 1 

This is the regime in which the term £5 gives the dominant contribution in Eq. (jlip. i.e. 

fi DE ~ 7/3r 2 . (17) 
Since r[ ~ 9ri/8 and r 2 ~ 3r 2 /8 from Eqs. (TT2"j) and (TTB"]) , the two variables T\ and r 2 evolve as 

7-1 oc a 9/s , r 2 cx a 3/8 . (18) 
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During the matter era we have 

w DE ~-l/8, w eS ~0. (19) 
If fir 2 > 0, then the scalar ghosts do not appear [IB, Ell- For the initial conditions with r 2 > we require that 

/8>0. (20) 
The conditions for the avoidance of Laplacian instabilities are automatically satisfied 0, [TtJ • 
• (ii) n = 1 

Along the tracker characterized by r\ = 1, there is a simple relation 

Ode = r 2 . (21) 

In the regime r 2 <C 1 one has r 2 ~ 6r 2 from Eq. (|13[) . so that the evolution of r 2 during the deep matter era is 
given by 

r 2 oc a 6 . (22) 
Along the tracker the evolution of the Hubble parameter can be analytically known as [4(| 



- ^O m0 (l + zf + Jl - Q m0 + (1+ / )6 (» m0 ) 2 , (23) 



H J 2 "' uv V 4 

where z — a®/ a— 1 is the redshift (ao is the scale factor today) , and fl m o is the density parameter of nonrelativistic 
matter today. Even in the presence of radiation and the cosmic curvature there exists an analytic form of H(z) 

Since p D E = 3M 6 /H 2 and p BE = -3M 6 (2 + w cS )/H 2 , it follows that 

2 2r 2 

W D E = -2 - Weft = — , W c ff = — ■ (24) 

r 2 + l r 2 + 1 

In the deep matter era (r 2 <C 1) we have i«de = — 2 and iy e g = 0, whereas at the dS solution wde = w e g = — 1. 
The conditions for the avoidance of ghosts and Laplacian instabilities have been derived in Refs. [lj| [l7[- The 
allowed parameter space in the (a,j3) plane is summarized in figure 1 of Ref. 

It should be pointed out that the tracker solution and the /^-dominant solution can be distinguished also by means 
of the following physical quantity. When c\ = one has dL/d<fi = 0, where L = ^—gC is the Lagrangian including the 
volume factor \J—g in Eq. (fT]). Then the conjugate momentum to the field <fi is conserved, so that dL/dcj) = constant. 
This constraint can be written as 

TV, 1 / 4 a 3 (1 - r x ) (2r 2 2 +3a ri 2 - 4/3 - 6 an + 5/3 n +5/3) , 
Prf, = t— j-. = constant . (25) 

On the tracker exactly vanishes 1 , whereas for the /^-dominant solution ~ 5/3 a 3 (r2/r} 1 ) 1 / 4 ^ 0. Although 
Eq. (|25|) is not independent of the others, we will use it later on in order to check the good convergence of numerical 
solutions. 

III. COSMOLOGICAL PERTURBATIONS 

Let us consider scalar metric perturbations ^, and x about the flat FLRW background [42|: 

As 2 = -(1 + 2*)di 2 - 2d i xdtdx i + a 2 (t)(l + 2<f>)5 ij dx i dx j , (26) 



1 The tracker can then be seen as the zero momentum solutions. 
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where we have chosen the spatial gauge such that the spatial metric is diagonal 2 . There is still a freedom to fix the 
temporal part £o of a vector associated with a scalar gauge transformation. Later on we will discuss this gauge degree 
of freedom in more detail. 

The field <j> is decomposed into the background and inhomogeneous parts, as <fi(t,x) = (j>(t) + 5(/>(t,x). In the 
following we omit the tilde for the background quantity. We shall use the action approach rather than taking the 
usual step of linearizing the equations of motion, as the former is more convenient from a computational point of view 
(indeed the scalar Lagrangian is quite compact compared to the tensorial equations of motion). Evidently the two 
approaches lead to equivalent equations of motion with the same set of solutions. For the perfect fluid, we cannot 
set w — in Eq. ([T]) identically from the beginning as its action would vanish. However, we can consider a fluid 
whose equation of state parameter w is different from zero, and only after we obtain the equations of motion for the 
perturbation variables, we will take the limit w — > 0. Since we are interested in the scalar perturbations of a cold 
fluid (temperature T = 0) with a barotropic equation of state of the form p m = wp m (in the limit w — > 0), the action 
approach of linear perturbation theory for a perfect fluid introduced in Refs. [43| can be simplified as follows. 

Let us define a 4- velocity of the perfect fluid, as u a = p^ 1 d a £. This relation defines £, which completely describes 
the scalar component of u a for the case under study. The normalization of the 4-velocity (u a u a = —1) implies that 
p = ' —g al3 da£dp£. Then at linear level we have 5fx = —5£ — pty. Now the quantity y/—gp m {^, s ) can be expanded 
at second order in the gravitational fields and in S£ and 5s. However, since the fluid is cold (T = 0), the first principle 
of thermodynamics, dp m = ndp — nTds, imposes that (dp m / ds)^, {d 2 p m / ds 2 )^, and (d 2 p m /dsdp) vanish 3 . In turn 
this implies that p m depends on p alone, and we do not need to expand the action in terms of 5s, the entropy 
perturbation field 4 . Since the action for matter is a function of p alone, as p m = p m (p), and 5p is a function of the 
perturbation variables 51 and ^ , the barotropic perfect fluid introduces only one new independent perturbed scalar 
field, 5£. Furthermore, since p m = np — p rn , it follows that dp m = pdn + nTds = pdn = ph5p/p at linear level. 
Using the field redefinition 5£ = —pv, where v is a velocity potential for the perfect fluid, after some algebra, we find 

v - 3Hwv = * + — — 5 , (27) 
1 + w 

where 5 = 5p rn /p m is the density contrast. 

In Fourier space the density contrast obeys the following equation of motion |44j 

6=-(l + w)(3$+£v-£^ , (28) 

where k is a comoving wave number. This equation follows from the first-order part of the continuity equation 
V^T^ = 0. In the following we derive other perturbation equations from the second-ord8er action for perturbations. 
In Appendix [3] we will show how to obtain Eq. (|28|) from the action approach. 

Let us consider nonrelativistic matter with w = 0. We introduce the gauge- invariant matter perturbation 5 m , as 

m 



5 + 3Hv . (29) 



Using Eqs. J28J and (J27J), it follows that 



■ k 2 ... 

5 m + 2HS m + —(y-x) = 3{Q + 2HQ) , (30) 
cr 

where Q = Hv — $. 

Let us now derive the coupled equations for metric/field/matter perturbations. Expansion of the action ([T]) at 
second order gives 

1(2) = 



h/=0Cp=L(*,*,X,ty,«). (31) 

Equation (|27|i shows that 5 and ^> are not independent fields. In the second-order action (|3T|) we have chosen the 
latter field. If we vary the Lagrangian (|3"Tj) with respect to "J, then we find a term p m (^f — v)/w that leads to an 



2 In the published PRD version there is an extra factor a(t) in front of the diX term. We corrected this typo. 

3 This also implies that n = n(/i). For example, fj, oc n w is enough to give p m = wp m . 

In the equations of motion approach, this result corresponds to the fact that, once Pm — "Wprni no entropy perturbations appear in the 
dynamical equations. 
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apparent singularity for w —> 0. However, Eq. ([27)) allows to reduce this term to —p m [5/(l + w) + 3Hv]. Taking the 
limit w — > at the end, we obtain the equation of motion 

k 2 k 2 (\ k 2 \ 

= A x * + A 2 5(j) - Pm v + A 3 — $ + + A 5 — x + -ciM 3 + A 6 — 5<j> - p m 6 = , (32) 

(X CL \ Zi (XI 

where the coefficients Ai are given in Appendix [B] For the theories up to the Lagrangian L 2 we have checked that, 
after substituting the relation v = "J [which comes from Eq. (|27p ] and using the background Friedmann equation for 
the coefficient of ^P, Eq. ([3"2"]) reduces to the usual 0-0 component of the perturbed Einstein equation obtained in the 
literature (see e.g. j45|). 

Variations of the second-order Lagrangian L with respect to $, X , and 8<j) result in the following equations 

k 2 ( k 2 3 
= + B 2 5(t> + + B A 5<i> + B 5 m + Bq—® + B 7 — + - Cl M 3 

a A \ a 1 2 

k 2 \ k 2 k 2 

B 8 —+B 9 )y + B 10 ^x + B u ^X + 3pmV = 0, (33) 

E x = Ci$ + C 2 5<j> + C 3 * + C A 6<j> + p m v = , (34) 

..... . k 2 

E H = Di$ + D 2 5<t> + L> 3 * + D A 84> + D 5 ^ + D Q — X 

k 2 \ k 2 ( k 2 \ k 2 

7 — + D 6 )& + D 9 —50 + D w — +Dn)V + D 12 — x = , (35) 



where the coefficients Bi , d , and Di are given in Appendix [Bj E$ corresponds to the spatial trace of the perturbed 
modified Einstein equations, E x to the momentum constraint, and Eg^ to the perturbed equation of motion for the 
Galileon field. Varying the Lagrangian L in terms of v gives rise to Eq. ([2"T]) . i.e. v = 'J for w = 0. 

Another constraint equation can be found by introducing the perturbation 7 in the metric (|26[) . as Sgij — 2a 2 didjj. 
Its equation of motion is given by 

E y = ^(5!$ + B 2 S<t> + B 3 <3> + Bjcj) + Bs* + B 9 ^>) + i Cl Af 3 (50 + Pm v = , (36) 

which does not contain 7 explicitly (as the quadratic terms in 7 in the action can be integrated out). Evaluating 
i?$ — 3i? 7 leads to the following equation 

£ T = B 6 $ + B 7 54> + B 8 * + B w X + B n X = . (37) 

The same relation can be derived by using the Bianchi identities, see Appendix lAl 

In what follows we choose the longitudinal gauge with x = 0. This choice, together with imposing 7 = 0, completely 
fixes the gauge. We also set c\ = to discuss the evolution of cosmological perturbations. 

IV. QUASISTATIC APPROXIMATION ON SUBHORIZON SCALES 

Let us consider the evolution of perturbations for the modes deep inside the Hubble radius. We derive the equations 
of matter perturbations and gravitational potentials under the quasistatic approximation on subhorizon scales (k ^> 
aH). This corresponds to the approximation under which the dominant contributions to the perturbation equations 
are those including k 2 /a 2 and S (46l. In the longitudinal gauge we obtain the following approximate equations 
from Eqs. ([32]), dSTJ) , and ([35]): 

k 2 

(A 3 $ + AeSt) ~ Pm 6 , (38) 

er 

B 6 <$> + B 7 5cf> + B s t> = 0, (39) 
L> 7 $ + D 9 5(f> + D 10 ^ ~ . (40) 



From Eqs. ([3"9"]) and (j4U)) one can express $ and 6(f> in terms of \E', i.e. 

A 3 D 9 - A 6 B 7 



9 ~ B 2 - B 6 D 9 * ' l41j 

■4 ~ A6 I S 6 ~ A f 7 ^, (42) 
B 7 — BqD 9 
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where we have used B% — A%, D-j = Bj, and Diq = Aq, Substituting Eqs. (|4Tj) and (|42l) into Eq. (|38|) . it follows that 

k 2 



* ~ -47rG cff p m (5 , (43) 



where 

2M 2 l (B% - B 6 D 9 ) 

Gcff = 2A 3 A 6 i? 7 -^ jD9 _ A 2 i?6 G' • (44) 

Since | j| is at most of the order of \HS\ in Eq. ([28]), one has |(fc 2 /a 2 )u| < |-ff<*|. This means that \Hv/S\ < {aH/kf < 1 
for subhorizon modes, which leads to S m ~ <5 in Eq. ()29|) . Since the r.h.s. of Eq. (|30|) can be neglected relative to 
the l.h.s., wc obtain the equation for the gauge-invariant matter perturbation for the modes deep inside the Hubble 
radius: 

5 m + 2H6 m - 4:TTG cS p m 5 m ~ , (45) 

where we have used the Poisson equation (|43j) . 

From Eqs. ((32) and (|4"3"| the field perturbation is given by 

3Af 2 (A 3 i? 7 -A 6 i? 6 ) ^qff^ 2 

* * UJ ^ ■ (46) 

We introduce the following quantity to describe the difference between the two gravitational potentials: 

$ A 3 D g -A 6 B 7 
7/= '^" BeZfr-B? ' (4?) 

where wc have used Eq. (|4"Tj) in the last approximate equality. We also define the effective gravitational potential 

$eff = (*-*)/2, (48) 

which is related with the deviation of the light rays in CMB and weak lensing observations [4l[. From the Poisson 
equation (|43|) it follows that 

^$ cff ~ -47rG cff ^-pmS ■ (49) 
Using the density parameter Q m = 8irGpm/(3H 2 ) and the relation S m — 8, we have 



3G off l + ?7 . f aH 
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$eS -"2"f 2 Um T 1 ' (50) 



The matter perturbation equation f|45[) can be written as 



2 + ^)s' m -^n m S m cO. (51) 



y m 2 g 

In the limit that </>—>■ and — > 0, the effective gravitational coupling G c s reduces to G. Hence, in the early cosmo- 
logical epoch, the General Relativistic behavior is recovered. If G e ff — G, then the evolution of matter perturbations 
during the deep matter era (H'/H ~ —3/2 and fl m ~ 1) is given by S m oc a cx i 2 / 3 . 

As the field velocity grows in time, G e g is subject to change. This leads to the modified growth rate of matter 
perturbations compared to the A-cold dark matter (ACDM) model. Unlike f(R) gravity [53] and Brans-Dicke theory 
[48| the effective gravitational coupling is independent of the wave number k. The quantity ?/ defined in Eq. (|47[) 
reduces to 1 for <j> — > and — > 0, but it deviates from 1 with the growth of <j> and cf>. In three different regimes 
characterized by (i) ri -C 1, r% "C 1, (ii) ri = 1, ri -C 1, and (iii) r\ = 1, r2 = 1, one can estimate G e ff and 77 as 
follows. 



• (i) n < 1, r 2 < 1 
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n = 1 + ( —P + —an ) r 2 + 0{r\) . (53) 



In this case we expand G e ff and rj about n — 0, r 2 — 0. Together with the use of the background equations (fT2)) 
and (fTU)) . we obtain 

= l+(^ + ^r 1 )r 2 + 0(rI) ! (52) 
V 8 M 16 

In this epoch the cosmological dynamics is dominated by the term £5. Since /3 > to avoid ghosts, we have 
that G c ff > G and 77 > 1 in this regime. Hence the growth rates of 6 m and ^eff are larger than those in the 
ACDM model. 

• (ii) n = 1, T2 -C 1 

Expanding G e ff and 77 about r 2 = 0, it follows that 

G cff i 291a 2 + 702/3 2 - 933a/3 + 20a - 84/3 + 4 ^. 2 , 

— = 1 + 2(10a-9/3 + 8) ^ + °^ ' (54) 

3(126a 2 + 306/3 2 - 405a/3 + 4a - 30/3) 2 

77 = 1 ; r 2 +0(r 9 ). (55 

' 2(10a-9^ + 8) - \ 2) V 1 

The evolution of G e ff and 77 depends on both a and /3. If a = 1.4 and (3 = 0.4, for example, we have 
Geff/G ~ 1 + 4.31r2 and 77 ~ 1 — 5.11r2, respectively. In this case G c g > G, but 77 is smaller than 1. 

• (hi) n = 1, r 2 = 1 

At the dS point we have 

Gee _ 1 /r^N 

G " 3(a-2/3)' 

r? = 1. (57) 

This means that 77 is not subject to change compared to the ACDM model. In Refs. [l|| [ljj it was shown 
that the viable parameter space consistent with the absence of ghosts and instabilities is confined in the region 
2/3 < a < 2/3 + 2/3. From Eq. ([56)) the effective gravitational coupling is constrained to be G e ff/G > 1/2. On 
the line a = 2/3, G e g goes to infinity (this includes the case in which only the terms up to £3 are present, i.e. 
a = P = 0). On the line a — 2(3 = 1/3, the effective gravitational coupling is equivalent to G. 

For the Lagrangian up to £3 one has r\ = 1, but the presence of scalar nonminimal couplings with curvature terms 
in £4 and £5 (such as R(V(f>) 2 /2) gives rise to the situation where 77 is different from 1. The above approximate 
formulas are useful to discuss the evolution of perturbations on the scales relevant to large-scale structures. 



V. NUMERICAL RESULTS 



In this section we present numerical results for the evolution of perturbations without employing the quasistatic 
approximation on subhorizon scales. The accuracy of the quasistatic approximation will be confirmed for the modes 
k > aH. 

Among the six equations (|32l - ((37|) . three of them are independent. For the numerical purpose we solve Eqs. (|32l) . 
([34]) . and (|37| together with Eqs. (|27j) and ([28]) . We have also confirmed that defined in Eq. (f25|) remains constant 
up to the accuracy of 10 -8 during the whole evolution. It is convenient to introduce the following dimensionless 
variables 

V = Hv, 5<p = 8(t>/(x dS M p {), (58) 

with A 1 = Ai/(HM$), A 2 = x dS A 2 /(HM pl ), A 3 = A 3 /M* h A\ = A 4 /(fl' 2 M p 2 1 ), A 6 = x dS A 6 /M ph B 6 = B 6 /M^, 
Bj = XdsB^/Mpi, and G4 = XdsCi/{H M p i). In the longitudinal gauge we obtain the following equations from 
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Figure 1: The equation of state iode versus the redshift z for the model parameters a — 1.37 and /3 = 0.44 with two different 
initial conditions: (A) n = 0.03, r 2 = 0.003, and (B) n = 0.999, r 2 = 7.0 x 10" 11 . The case (B) corresponds to the tracker 
solution, whereas in the case (A) the solution approaches the tracker around today. The present epoch (z = 0) is shown as a 
dotted line (we also draw the dotted line in other figures). 



Eqs. H32J, (jHJ), dHTJ), (EHD, and dHJ): 

* = -(-B 6 $ + B 7 5p>)/A 3 , (59) 

$' = [(3i 4 i 6 B6 + A^Be - 3A 2 A 6 k 2 /{aH) 2 - 9A 6 B 6 n m )$ 

+ {3A 2 A 3 C A + 3A 4 A 6 B 7 + A X A 2 B 7 - 9A 6 B 7 Q m - 3A 3 A 2 k 2 / \aH) 2 )Sp + 9A 3 A 6 il m S + 9A 2 A 3 Q m V] 
x[3A 3 (A 1 A 6 - A2A3)}- 1 , (60) 

V = -[{A\B 6 + 3A 3 A 4 B 6 - 3A 3 3 k 2 /{aH) 2 - 9l 3 B 6 O m )$ 

+{A\B 7 + 3A 3 AiB 7 + 3AxA 3 Ci - 3A 2 3 A 6 k 2 /(aH) 2 - 9A 3 B 7 Q, m )5<p + 9A 2 3 a m S + 9AxA 3 9, m V] 

x[3l 3 (Iil 6 - i 2 i 3 )r 1 , (61) 
5' = -3t>' -k 2 /{aH) 2 V , (62) 
V' = (H'/H)V + ^ , (63) 

where we have used B% = A 3 , C\ = A 3 , C 2 = A§, and C 3 = — Ai/3. The time-dependent coefficients Ax, Bq, c.t.c. 
can be expressed by using the variables a, /3, fx, and r 2 . Solving the perturbation equations (|59l) -( l63j) together with 
the background equations (jT^j) and P^|) . we find the evolution of VP, $, S(p, 5, and V. 

During the deep matter era in which the field perturbation is negligibly small, the evolution of ^, <£>, <5, and V 
is similar to that in GR. The initial conditions are chosen to satisfy <£>' = and S(p' = 0. From Eqs. and (|6"Tj) 
this gives two constraints on the variables <£>, V, 6ip, and S. For given Sip and S, the initial conditions of $ and V 
are determined accordingly. For the modes deep inside the Hubble radius, Eq. PS]) provides a relation between Sip 
and (5. We choose the initial condition S = 10 -5 , but this can be chosen to be any value (as long as the amplitudes 
of perturbations do not matter). We note that the relation pp]) cannot be used for the wavelengths larger than the 
Hubble radius. Apart from the modes deep inside the Hubble radius at the onset of integration, we adopt the initial 
condition Sip = 0. Provided that Sip is small initially, the dynamics of perturbations is similar to that in the case 
Sp = 0. For large initial values of Sip, the field perturbation tends to oscillate for small-scales modes (associated with 
large wave numbers k). This situation is similar to what happens in the generalized Galileon model [20(. 

We identify the present epoch (the redshift z = 0) to be J7de = 0.72. In order to find the evolution of the quantity 
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Figure 2: Evolution of the perturbations for a = 1.37 and /3 = 0.44 with the background initial conditions n = 0.03 and 
T2 = 0.003 (corresponding to the case (A) of Fig. [T]). (Left) 5 m /a versus z for the wave numbers (a) k = 300etoiio, (b) 
k = 30ao-Ho, and (c) k = 5ao-Ho. (Right) <J> e ff versus z for the wave numbers (a) k = 300ao-ffo, (b) = lOao-ffo, and (c) 
k — 5aoHo- Note that Sm/a and <fr cf f are divided by their initial amplitudes Sm{U) / a{U) and $ e ff(ti), respectively, so that their 
initial values are normalized to be 1. The bold dotted lines show the results obtained under the quasistatic approximation on 
subhorizon scales. The choice of initial conditions for perturbations is explained in the text. 

k/(aH), we also integrate the following equation 

b' = b(l + H'/H) , (64) 

where b = aH. The scales relevant to the linear regime of the galaxy power spectrum are 30 < k/(aoHo) ;$ 600, where 
the subscript "0" represents present values. Note that the above upper limit corresponds to k = 0.2/iMpc -1 , where 
h = 0.72 ± 0.08 |4{|. For the scales k/(ctoHo) > 600 the nonlinear effect becomes crucially important. Meanwhile, the 
wave numbers relevant to the ISW effect in CMB anisotropies correspond to the large-scale modes with k/(aoHo) = 
0(1). 

In Ref. [4(| two of the present authors placed observational constraints on the Galileon model ([1]) from the back- 
ground expansion history of the Universe. The combined data analysis of the type la supernovae (Constitution and 
Union2 sets), the CMB shift parameters (WMAP7), and the baryon acoustic oscillations (BAO) show that the tracker 
solution described by Eq. (|23l) is not favored, but the solutions that approach the tracker at late times can be compat- 
ible with the background observational constraints. As illustrated in Fig. [T] the equation of state u>de for the tracker 
[case (B)] evolves from —2 (deep matter era) to —1 (dS epoch). The solution (A) in Fig. [1] starts to evolve from the 
value wde = — 1/8 and then it enters the tracking regime around the present epoch. In the flat FLRW background 
the best-fit model parameters are a = 1.411 ± 0.056, /3 = 0.422 ± 0.022 (Constitution + CMB + BAO, 68% C.L.), 
and a = 1.404 ± 0.057, /? = 0.419 ± 0.023 (Union2 + CMB + BAO, 68 % C.L.). 

A. Case of late-time tracking solutions 

Let us consider the evolution of perturbations for the late-time tracking solutions. In Fig. [5]we plot S m /a and $ c ff 
versus the redshift z with the model parameters and initial conditions corresponding to the case (A) of Fig. [1] For 
the mode k = 300ao-ffo we nn d that the full numerical result shows excellent agreement with that obtained under the 
quasistatic approximation on subhorizon scales. The difference starts to appear for the modes fc/(ao-ffo) < O(10). 
The left panel of Fig. [2] shows that, on larger scales, the growth of S m tends to be less significant. For the modes 
k ^> ao-ffo the matter perturbation evolves faster than a during the matter era (i.e. faster than in the case of GR, 
5 m oc a). The growth of S m /a turns into decrease after the Universe enters the epoch of cosmic acceleration. 
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Figure 3: Variation of the growth index 7 of matter perturbations for a = 1.37 and /3 = 0.44 with the mode k — 300ao-Ho. The 
initial conditions of n and 7-2 for the cases (A) and (B) are the same as those given in the caption of Fig. [T] 



The case (A) in Fig. [^corresponds to the same model parameters and initial conditions as those given in the numerical 
simulations of Fig. [2] Unlike the ACDM model in which 7 is nearly constant (~ 0.55 [5l|) in the low-redshift regime 
< z < 1, the variation of 7 is significant for the solutions that approach the tracker at late times. Moreover the 
growth index today for the mode k = 300ao-ffo is 7o = 0.35, which is quite different from that in the ACDM. For the 
wave lengths relevant to large-scale structures we find that the growth index today exhibits almost no dispersion with 
respect to k. This property is different from viable f(R) dark energy models in which 70 can be dispersed [52j . If we 
choose smaller initial values of r\ (i.e. the later tracking), we find that 70 tends to be smaller, e.g. 70 = 0.23 for the 
initial conditions r\ = 0.01, r 2 = 0.003 with the model parameters a = 1.37, f3 = 0.44. 5 

The right panel of Fig. [2] shows that, unlike the ACDM model, the effective gravitational potential $ c ff changes in 
time even during the matter era for the modes k 3> aH. This can be understood as follows. For the initial conditions 
corresponding to the late-time tracking solutions the effective gravitational coupling G c ff and the anisotropic parameter 
j] are given by Eqs. ([5"2"j) and ([5"3")l . respectively, before reaching the tracker. Since G e g/G ~ 1 + 255/3r 2 /8 > 1 and 
i] ~ 1 + 129/3^/8 > 1 in this regime, we have the larger growth rates of $ e ff and S m relative to those in GR. In particular 
the term (G c ff/G)(l + 77) /2 in Eq. ((50)) is larger than 1, which leads to the additional enhancement of $ e ff to the 
growth coming from S m . Note that in f(R) gravity [13] and in Brans-Dicke theory [48| one has (G e ff/G)(l + ??)/2 = 1, 
so that the evolution of 5 m is directly related with the variation of $ c ff ■ In Galileon gravity the unusual behavior of 
the anisotropic parameter rj leads to the nontrivial evolution of perturbations. For the model parameters a = 1.37 
and (3 = 0.44, Eq. (|56| gives G c s — 0.68G at the dS fixed point. As we see in Fig. [2l $ e ff begins to decrease at some 
point after the matter era. 

For the large-scale modes relevant to the ISW effect in CMB anisotropies, i.e. k/(a^H^) < 10, the effective 
gravitational potential is nearly constant in the early matter-dominated epoch, see the cases (b) and (c) in the right 
panel of Fig. [2] However, $ e ff exhibits temporal growth during the transition from the matter era to the epoch of 
cosmic acceleration. Note that in the ACDM model $ e ff decays without the temporal growth after the matter era. 
The characteristic variation of $ e g in the Galileon model should leave observational signatures in CMB anisotropies 
as the ISW effect. 




(65) 



For the dark energy models with constant u>de and 7. the current observations still allow the large parameter space of 7 ranging from 
0.2 to 0.6. [H. 
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Figure 4: Evolution of (i) 5 m /a and (ii) $ c ff versus z for a = 1.37 and f3 = 0.44 with the initial conditions r\ = 0.999 and 
T2 = 7.0 x 10" 11 [corresponding to the case (B) of Fig. [T]. Both <5 m /a and 5> e ff are divided by their initial amplitudes. In this 
case the background cosmological solution is on the tracker from the onset of integration. The solid line shows the evolution of 
perturbations for the mode k — 5aoHo, whereas the bold dotted line and bold dashed line represent the result derived under 
the quasistatic approximation on subhorizon scales. 

The above numerical results correspond to the fixed values of a and /3 (a — 1.37, j3 — 0.44). If we use the bounds 
a = 1.404±0.057, (3 = 0.419±0.023 constrained by Union2+CMB+BAO data sets, the effective gravitational coupling 
at the dS point is restricted in the range 0.5G < G c s < 0.72G. The bounds coming from Constitution+CMB+BAO 
data give the similar constraint, i.e. 0.5G < G s < 0.71G. For the model parameters close to the upper limit 
a = 2(3 + 2/3 of the allowed parameter space (i.e. G e s is close to 0.5G at the dS point), the parameter r] tends to 
show a divergence during the transition from the matter era to the dS epoch. If G ff is larger than 0.66G, we find 
that such divergent behavior can be typically avoided. 

For the parameters a and (3 constrained observationally, the values of 70 are usually less than 0.4 for the late-time 
tracking with the minimum values of wde larger than —1.3. In addition the later tracking leads to smaller values 
of 70 with larger variations of 7. This property of the Galileon model can be clearly distinguished from the ACDM 
model. 



B. Case of early-time tracking solutions 

If the solutions arc already close to the tracker in the early matter-dominated epoch, the evolution of perturbations 
is different from that discussed above. For the tracker one has r2 oc a e during the matter dominance, which is much 
faster than the evolution T2 oc a 3 / 8 in the regime ri -C 1 and r2 < 1. This means that T2 is very much smaller than 
1 for the redshift z > 1. Then the second terms on the r.h.s. of Eqs. (ji>4"| and (|55|) are suppressed relative to the first 
terms until recently. In Fig.@]we plot the variation of 5 m /a and $ e ff for ol = 1.37 and (3 = 0.44 with the background 
initial conditions corresponding to the case (B) in Fig. [TJ The solid line shows the evolution of perturbations for the 
mode k = §o,qHq, whereas the bold dotted line and bold dashed line represent the result derived under the quasistatic 
approximation on subhorizon scales. In this case the evolution of the large-scale mode with k = ba^Ho is similar to 
that for the modes deep inside the Hubble radius. 

Figure [4] shows that the growth of matter perturbations in the deep matter era is almost identical to that in the 
ACDM model. In the regime r\ ~ 1 and r2 < 1 one has G e s/G ~ 1 + 3. 2^2 and 77 ~ 1 — 3.7r2, both of which are 
very close to 1. Recall that at the dS point G e s/G ~ 0.68 and i] — 1. Since 77 is very close to 1 and G e g is different 
from G only around the dS solution, the evolution of perturbations is milder than that corresponding to the late-time 
tracking solutions. The case (B) of Fig. [3] shows 7 versus z for the mode k = 300ao-ffo with the initial condition 
corresponding to the tracker solution. In this case the growth index today is found to be 70 — 0.47 with the variation 
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Figure 5: Evolution of perturbations for a — and P = with the initial conditions ri = 0.05 and = 0.001. The lines 
correspond as follows: (Al) 8 m /a for the mode k = 300ao-Ho, (A2) <J? e ff for the mode k — 300ao-Ho, (Bl) 5 m /a for the mode 
k — 5aoHo, and (B2) $ G fr for the mode k = 5aoHo- Both <5 m /a and <E> e ff are divided by their initial amplitudes. The bold 
dotted line and bold dashed line correspond to the results obtained under the quasistatic approximation on subhorizon scales. 

of 7 in the low-redshift regime. In Fig. 0] we also find the small temporal growth of $ c ff from the end of the matter 
era to the epoch of cosmic acceleration. This can also give rise to a nonncgligible contribution to the ISW effect in 
CMB anisotropics. We recall, however, that the tracker solution is not favored at the background level by the joint 
analysis of observational data. 



Let us finally consider the model in which only the terms up to £3 arc present, i.e. a = /3 = 0. Although such a 
model is disfavored observationally, it is of interest to study what happens for the evolution of perturbations due to 
the peculiar divergence of G e ff at the dS point (which occurs along the line a = 2/3). 

First of all we have numerically found that the late-time tracking behavior as in the case (A) of Fig. [T]is not easy 
to be realized even by choosing many different initial conditions. This is the main reason why the combined data 
analyses at the background level do not favor such model parameters. For a = (3 = one has 77 = 1 independent of 
the values of r\ and r2, whereas G s — G for r\ -C 1, r<z "C 1 and G g ~ 1 + r^/A for r\ = 1, r-z -C 1. 

In Fig. [5] we plot the variation of S m /a and $ e ff with the initial conditions n = 0.05 and r2 = 0.001 for two different 
wave numbers: k = 300ao-£/o an d k = 5cto-ffo- I n this simulation the solution reaches the tracker region (ri ~ 1) 
around the redshift z = 5, e.g. n — 0.99 and r2 = 4.14 x 10 -4 at z = 5.34. Hence the evolution of perturbations for 
z > 5 is similar to that in the ACDM model. The effective gravitational coupling G e g increases with the growth of ?'2 
for z < 5. In this case the numerical value of G e ff today is G e ff — 1.94G with r-i — 0.72. For the mode k = 300ao^o 
the growth index today is found to be 70 — 0.43, which is smaller than the value in the ACDM model. 

In Fig. [S] the effective gravitational potential for the mode k = 300ao-ffo grows around the present epoch in spite 
of the decrease of S m /a. This can be understood by the fact that the term (G e ff/G)(l + r])/2, which is equivalent to 
Geff/G for a = ft = 0, continues to grow toward the dS solution. While the decrease of the density parameter fi m 
overwhelms the increase of G e s/G in Eq. (|5"Tj) for matter perturbations, the presence of the additional term (aH/k) 2 
in Eq. (|50[) stimulates the growth of $ c ff around today. Numerically we find that the growth of $ c ff eventually turns 
into the decrease in future as 6 m decays sufficiently (which is not shown in Fig. [5]). For k = 5cioHq the growth of $ c ff 
around today is milder relative to the modes deep inside the Hubble radius. 

The anticorrelation between 5 rn /a and $ c ff seen in Fig. [5] is similar to the one found in the generalized Galileon 
model |2(j. This anticorrelation mainly comes from the fact that G e g/G continues to grow toward the dS point and 
that the parameter 77 does not become smaller than 1 to compensate the growth of G e g/G (unlike /(i?) gravity and 
Brans-Dicke theory) . In the numerical simulations of Fig. [2] the effective gravitational coupling is larger than G in 



C. Case of a = /? = 
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the regime r\ <C 1 and T2 -C 1, but it finally approaches the asymptotic value G e ff — 0.68G. In this case both S m /a 
and $ e ff start to decay by today, whose property is different from that in the case a = f3 = 0. Since the observational 
constraints on a and /3 do not allow the region close to the line a = 2/3, the anticorrelation between S m /a and $ e ff 
does not typically occur for the viable model parameters. 

VI. CONCLUSIONS 

We have studied the dynamics of cosmological perturbations in the Galileon model whose Lagrangian satisfies 
the Galileon symmetry <9 M — ¥ d^4> + in the flat space-time. In this theory there exists stable dS solutions with 
</> = constant. In the deep matter era the General Relativistic behavior can be recovered because of a small field 
velocity, but the Universe finally enters the epoch of cosmic acceleration after the growth of cf>. 

Before reaching the dS attractor, the solutions approach a tracker characterized by r\ = 1 (i.e. the evolution </> oc 
H^ 1 ). The tracking epoch depends on the initial conditions of the background cosmological variables. The combined 
data analysis based on the background expansion history of the Universe (SN Ia+CMB shift parameters+BAO) favor 
the late-time tracking around today. This corresponds to the initial conditions r\ <C! 1 and <C 1 in the deep matter 
era, so that r\ approaches 1 in the low-redshift regime (z < 1). 

The study about the evolution of cosmological perturbations can allow us to distinguish the Galileon model from the 
ACDM model further. In the presence of a nonrelativistic perfect fluid, we have derived full perturbation equations 
for the model described by the action ([1]). In spite of their complexities it is possible to obtain simpler equations 
for the matter perturbation 6 m and the effective gravitational potential $ c ff under a quasistatic approximation on 
sub horizon scales. 

The two important quantities associated with the growth of S m and $ c ff are the effective gravitational coupling G ff 
and the anisotropic parameter rj = — $/^. In three distinct regimes we have approximately derived the expressions 
of G e ff and r], see Eqs. (|52 p -(|57 p . For the initial conditions that lead to the late-time tracking behavior, the solutions 
start from the regime n <C 1, ^2 <C 1 with positive /3 (which is required to avoid ghosts). In this regime Eqs. (|52[) 
and (JS3J) show that G c g/G ~ 1 + 255/3^/8 > 1 and r) 1 + 129/3^/8 > 1. This gives rise to the larger growth rates 
of 8 m and $ e ff relative to those in the ACDM model. 

In f(R) gravity and Brans-Dicke theory the anisotropic parameter r\ is less than 1, so that the combination 
(G e ff/G)(l + 77) /2 remains to be 1. In the Galileon model the fact that rj > 1 in the regime v\ <C 1, r2 <C 1 
leads to the further growth of $ e ff for subhorizon modes. At the dS fixed point we have that G e ff/G = l/[3(a — 2/3)] 
and 77 = 1. For the parameters of a and (3 observationally constrained at the background level (40[,tne effective 
gravitational coupling G e ff is smaller than 0.72G. Then the perturbations S m /a and $ e ff turn into decrease around 
today, as we see in the numerical simulation of Fig. [5J 

For the solutions that approach the tracker at the early epoch of the matter era, G e ff and 77 are approximately given 
by Eqs. (|54[) and (|55[) . respectively, in the regime n = l,f2 <C 1. Since the evolution of r2 in the tracking regime is 
very fast (r2 oc a 6 ), the correction term /3t~2 to G e ff and rj becomes important around the present epoch in which 
grows to the order of 1. Hence the deviation from the ACDM model for the evolution of S m and $> c s appears only at 
late times, as we see in Fig. 2J The early tracking behavior is however disfavored from joint observational constraints 
at the background level. 

For the late-time tracking solutions we have found that the growth index 7 of matter perturbations rapidly changes 
in the low-redshift regime [as in the case (A) of Fig. [3] . The values of 7 today on the scales relevant to large-scale 
structures are typically smaller than 0.4 for the model parameters constrained by the observations at the background 
level. As the tracking occurs at a later epoch, the presents values of 7 tend to be smaller. On large scales relevant 
to the ISW effect in CMB anisotropics the effective gravitational potential $ G ff exhibits a temporal growth around 
the present epoch. These properties can potentially restrict the allowed region constrained by the background cosmic 
expansion history further. It will be of interest to constrain the Galileon model from the combined data analysis of 
large-scale structures, CMB, and weak lensing. 
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Appendix A: Bianchi identities 

We confirm the consistency of perturbation equations by using the Bianchi identities. We write the action (fTJ) in 
the form 



<Txy/=g£, (Al) 
where C is the total Lagrangian. Variation of the action (|A1|) with respect to g a p gives 

5S = J tfxy/^jjZ^Sgcf). (A2) 

Then the Bianchi identities lead to 

V /3 £ Q ' 3 = d^ 13 + r^s^ + r^E QA = , (A3) 

where T"^ is the Christoffel symbol. 

We perturb Eq. (|A3|) at first order. Then we find 



dpSE* + r^<5I] A/3 + rJj<J£ a * = , (A4) 
where we have used the fact that the background equation satisfies S Q/3 = 0. For the case a = we have 

d 6X 00 + d l 8Yp l + T%SY xp + r^(5S 00 + rf fl JE w = . (A5) 



00 — u > x Oi ~ U J x ij — " 1JU >) ' Tq = 0, Tq^- = H8ij 7 T* k 

d 5Z oa + d l 5Y m + a 2 H 5 ^ U + = . (A6) 



Using the relations r° = 0, rg j = 0, r°- = a 2 H8 ih = 0, TL = HSij, T) k = 0, Eq. reduces to 



In our Galilcon model the condition (|A6|) leads to the following relation 

Si = Ess, + 3HE^ - HE<s, --^E x = 0, (A7) 

where Ey corresponds to Eq. ([52]). to Eq. ([55]). and £ x to Eq. ([51]). Combining Eq. ([35]) with Eq. ([AT]) and 
inserting the background equations of motion, we find 

(k 2 k 2 A 

E 1 - <j>E S 4 = p m ^x-3$-^-«-c$ =0, (A8) 



which matches with Eq. ([28]) for «; = 0. 

Equation ([37]) can be also derived from the spatial part of the Bianchi identities. It follows that 



2 



E x + 3HE X - ^E* = I k —E^ = , (A9) 



where E y corresponds to Eq. ([57]). 



Appendix B: Coefficients of perturbation equations 

In this Appendix we show the coefficients of the perturbation equations ([32]) . ([33]) . (|34|) . and ([35l) . 

Eg. m 

A x = 6HAI 2 ! - 3c 3 <p 3 /M 3 + 45c 4 H^/M 6 - 63c 5 H 2 <j> 5 /M 9 , (Bl) 

A 2 = c 2 - 9c 3 H^) 2 /M 3 + 90c4,H 2 (j) 3 /M 6 - 105c 5 £T 3 <p 4 / M 9 , (B2) 

A 3 = 27\/ 2 ! + 3c 4 cj) 4 /M 6 - 6c 5 H<j) 5 /M 9 , (B3) 
A A = 2p m - 9H 2 M^ - Cl M 3 (j)/2 - 3c 2 2 /2 + 15c 3 i7</> 3 /A/ 3 - 315c 4 iJ 2 4 /(2M 6 ) + 189c 5 H 3 <j> 5 /M 9 , ( B4 ) 

A 5 = -2HM^ + c 3< /) 3 /M 3 - WciH^/M 6 + 2lc 5 H 2 <j) 5 /M 9 , (B5) 

-46 = -c 3 4> 2 /M 3 + 12 Ci H<j) 3 /M 6 - 15c 5 H 2 ^/M 9 . (B6) 
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Eg. gj 

B 1 = 3A 3 , (B7) 

B 2 = 3A 6 , (B8) 

B 3 = 18HM^ + 27c 4 B</> 4 /Af 6 + 36c 4 </> 3 0/A/ 6 - 54c 5 H 2 <j) 5 /AI 9 - 90c 5 H<p 4 4>/AI 9 - 18c 5 B</> 5 /Af 9 , (B9) 
B 4 = -3c 2 <j) - 6c 3 00/M 3 + 36c 4 B0 3 /A/ 6 + 108c 4 B <j) 2 4> / 'AI 6 + 54c 4 B 2 3 /Af 6 

-90c 5 H 3 ^/M 9 - 180c 5 H 2 <j) 3 4>/M 9 - 90c 5 HH^/AI 9 , (BIO) 

B 5 = -A u (Bll) 

B 6 = 2M^ - c^/AI 6 - 6c 5 <j) 4 4>/M 9 , (B12) 

B 7 = 12c 4 <p 2 4>/M 6 + 4c 4 B</> 3 /A/ 6 - 6c 5 H<p 4 /M 9 - 6c 5 B 2 4 /A/ 9 - 24c 5 B> 3 ^/A/ 9 , (B13) 

B 8 = A 3 , (B14) 
B 9 = -9H 2 M 2 l - 6iJA/p! - 3p TO + 3ciA/ 3 < ? !)/2 + 3c 2 2 / 2 + 9c 3 </. 2 0/A/ 3 - 135c 4 B 2 07(2A-f 6 ) 

-18Oc 4 B0 3 0/A/ 6 - 45c 4 B</> 4 /A/ 6 + 315c 5 H 2 ^4>/AI 9 + 126c 5 H 3 <j) 5 / M 9 + 126c 5 HH<j) 5 /M 9 , (B15) 

Bio = -A 3 , (B16) 

Bii = -2BAf 2 i - 3c 4 B</> 4 /A/ 6 - 12c 4 3 <?i/A/ 6 + 6c 5 H<p 5 /AI 9 + 3Oc 5 B0 4 0/A/ 9 + 6c 5 H 2 <j) 5 /M 9 . (B17) 

Eg. AMD 

d=A 3 , C 2 =A 6 , C 3 = -A 1 /3, (B18) 

C 4 = -C2</> + 3c 3 B0 2 /M 3 - m Ci H 2 (j) 3 /M 6 + lbc 5 H 3 <j) 4 /M 9 . (B19) 

Eg. (EH) 

£>i = SA e , (B20) 

B 2 = c 2 - 6c 3 H<j)/M 3 + M Ci H 2 <i) 2 /M 6 - 60c 5 H 3 4> 3 /M 9 , (B21) 
D 3 = 3c 2 (j) - 18c 3 H<j) 2 /M 3 - 6c 3 (j)4>/M 3 + lO8c 4 B0 2 0/A/ 6 + 162c 4 B 2 c/> 3 /A/ 6 + 36c 4 H<j) 3 /M 6 

~90c 5 (HH(j)' l /M 9 + 2H 3 <j) 4 /M 9 + 2H 2 $ 3 $/AI 9 ) , (B22) 
B 4 = 3c 2 H - 6c 3 B</>/M 3 - 18c 3 H 2 <j)/AI 3 - 6c 3 B0/AB 3 + 162 Ci H 3 <j> 2 / AI 6 + lO8c 4 BB0 2 /Af 6 

+lO8c 4 B 2 00/M 6 - l8Qc b {H 3 4> 2 4>/M 9 + H 2 H<p 3 /M 9 + iJ 4 ^> 3 /M 9 ) , (B23) 

D 5 = -A 2 , (B24) 

B 6 = -A 6 , (B25) 

B 7 = B 7 , (B26) 
B 8 = 3c 1 M 3 /2 + 3c 2 (/> + 9c 2 H4> - 27c 3 H 2 <j) 2 /AI 3 - 18c 3 H <j)4> / M 3 - 9c 3 H<j> 2 /AI 3 + l62c 4 H 2 <j) 2 j> / M 6 

+162C4,H 3 4> 3 /M 6 + lQ8c 4 HH<p 3 /M 6 - \80c 5 H 3 j) 3 4> / M 9 - 135c 5 (B 2 H ^ / 'AI 9 + B 4 4 /Af 9 ) , (B27) 
D g = c 2 - 4c 3 H<j)/M 3 - 2c 3 (j>/M 3 + 26c 4 B 2 </> 2 /AI 6 + 12c 4 B</> 2 /A/ 6 + 24c 4 B<^/Af 6 

-36c 5 H 2 <p 2 4>/Al 9 - 24c 5 (H 3 <j) 3 /M 9 + HH<j) 3 /AI 9 ) , (B28) 

Bio = A 6 , (B29) 
Bn = ciM 3 /2 - 3c 2 H(j) - c 2 (j) + 27 c 3 H 2 <j) 2 /AI 3 + 18c 3 H^4>/AI 3 + 9c 3 H(j) 2 /M 3 - 18Oc 4 BB0 3 /M 6 

-27Oc 4 (B 2 2 0/A/ 6 + H 3 <p 3 /M 6 ) + 420c 5 B 3 (/) 3 ^/Af 9 + 315c 5 (B 4 (/) 4 /Af 9 + H 2 H4> A /M 9 ) , (B30) 
B12 = -c 2 </> + Ac 3 H4> 2 /M 3 + 2c 3 <j)j)/AI 3 - 30c 4 B 2 ^ 3 /Af 6 - 36c 4 B0 2 0/A/ 6 - 12c 4 B0 3 /A/ 6 

+3Oc 5 (BB0 4 /7l/ 9 + H 3 4> 4 /AI 9 + 2H 2 4> 3 '<j)/Ad 9 ) . (B31) 
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